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We introduce the Korringa-Kohn-Rostocker non-local coherent potential approximation (KKR- 
NLCPA) for describing the electronic structure of disordered systems. The KKR-NLCPA systemat- 
ically provides a hierarchy of improvements upon the widely used KKR-CPA approach and includes 
non-local correlations in the disorder configurations by means of a self-consistently embedded clus- 
ter. The KKR-NLCPA method satisfies all of the requirements for a successful cluster generalization 
of the KKR-CPA; it remains fully causal, becomes exact in the limit of large cluster sizes, reduces 
to the KKR-CPA for a single-site cluster, is straightforward to implement numerically, and enables 
the effects of short-range order upon the electronic structure to be investigated. In particular, it is 
suitable for combination with electronic density functional theory to give an ab-initio description 
of disordered systems. Future applications to charge correlation and lattice displacement effects 
in alloys and spin fluctuations in magnets amongst others are very promising. We illustrate the 
method by application to a simple one-dimensional model. 

PACS numbers: 71.15.Ap, 71.23.-k, 71.20.Be, 71.15.mb 



I. INTRODUCTION 

Over the past 30 years or so the coherent potential ap- 
proximation-* (CPA) has proved to be a generally reliable 
method for dealing with disordered systemsi^ However, 
being in essence a single-site mean-field theory^ the CPA 
fails to take into account the effect of non-local potential 
correlations due to the disorder in the environment of 
each site and hence leaves much important physics out 
of consideration. Consequently, considerable effort has 
been spent in trying to find a way of improving it system- 
atically by a multi-site or cluster generalization. Surpris- 
ingly this has turned out to a very difficult problem^ and 
a viable solution has been proposed only recently. The 
new method has emerged from the Dynamical Cluster 
Approximation 6 - 7 - 8 (DCA) which was directed originally 
at describing dynamical spin and charge fluctuations in 
simple Hubbard models of strongly-correlated electron 
systems. Recently its static limit has been adapted by 
Jarrell and Krishnamurthy for a simple tight-binding 
model of electrons moving in a disordered potential^ The 
same problem was investigated by Moradian et al^fi In 
this paper we develop the ideas behind this approxima- 
tion further and demonstrate how they can be combined 
with realistic, ab-initio descriptions of systems of inter- 
acting electrons in disordered systems. 

Because the language of our multiple scattering the- 
ory is so different from that of the context in which the 
DCA is usually deployed, we elaborate on this relation- 
ship. Firstly, we note that the DCA was invented to 
describe short-range correlations within the framework 
of the Dynamical Mean Field Theorjiii (DMFT) of spin 
and charge fluctuations in many-electron systems. Sec- 
ondly, we recall that that the DMFT can be regarded as 
the dynamical generalization of the CPA for the Hubbard 
'alloy analogy' problem liSiiii Thus, in the light of these 



remarks, it is natural to investigate the static version of 
the DCA as a generalization of the CPA which includes a 
description of short-range order. Indeed, Jarrell and Kr- 
ishnamurthjai already studied the problem of electronic 
structure in random alloys from this point of view. In 
this paper we tackle the same problem using an identical 
conceptual framework but a very different description of 
the electrons afforded by multiple scattering theoryii2ii£i 
To make the above remarks more specific we would like 
to highlight two of the principal differences between our 
treatment of the problem and that of Ref. 0. Firstly, 
we do not make use of a tight-binding model Hamilto- 
nian but solve, numerically, a Schrddinger equation in 
each unit cell and match the 'out-going wave' solution to 
the incoming waves from all the other unit cells. This 
is known as the multiple scattering approach-l^i^ to the 
problem of electronic structure in solids and is the foun- 
dation of the Korringa-Kohn and Rostoker (KKR) band 
theory method. Consequently, the principal virtue of our 
formalism, as opposed to those based on tight-binding 
model Hamiltonians is that it prepares the ground for 
first-principles calculations based on density-functional 
theories i& The second difference is a formal consequence 
of the first. In multiple scattering theories the object of 
interest is not the self-energy and the diagrammatic lan- 
guage of perturbation theory is not used. We will show 
that the quantities that play the role of the self-energy 
in multiple scattering theory are the effective scattering 
amplitudes t and effective structure constants G(R y ) 
which are also the natural concepts in effective medium 
theories. 17 In short, these formal reasons fully account for 
the fact that we do not base our arguments on 'restoring 
momentum conservation' and introducing approximate 
Laue functions to renormalize diagrams but construct 
our theory in terms of real and reciprocal space clus- 
ters. Nevertheless, we believe that our final algorithm 
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described in Sec. Ill El is equivalent to those investigated 
by Jarrell and Krishnamurthy. 9 Our aim in reformulating 
the problem is to facilitate the deployment of the method 
as a first-principles calculation, in other words to develop 
a non-local KKR-CPAi& 

In brief, our KKR-NLCPA method introduces new ef- 
fective structure constants and this enables us to de- 
fine an effective medium which includes non-local poten- 
tial correlations over all length scales. Using a 'coarse- 
graining' procedure inspired by the DCA we can then de- 
rive a self-consistent 'cluster generalization' of the KKR- 
CPAi& which determines an approximation to this effec- 
tive medium by including non-local correlations up to the 
range of the cluster size. The KKR-NLCPA satisfies all 
of the requirements for a successful cluster generalization 
of the KKR-CPA as listed by Gonis£ In particular, the 
KKR-NLCPA becomes exact in the limit of large clus- 
ter sizes where it includes non-local correlations over all 
length scales, and recovers the KKR-CPA for a single- 
site cluster. The method is fully causal, allows the effects 
of short-range order to be modelled, and can be imple- 
mented numerically for realistic systems. 

The outline of this paper is as follows. In the next sec- 
tion we describe the formalism for the KKR-NLCPA. We 
explain our KKR-NLCPA algorithm and show how to in- 
clude short-range order. We describe in more detail how 
to carry out the coarse-graining with reference to simple 
cubic, body-centered cubic and face-centered cubic lat- 
tices. Finally we explain how to use the KKR-NLCPA 
formalism to calculate observable quantities such as the 
configurationally-averaged density of states in prepara- 
tion for DFT calculations. In order to illustrate the im- 
provements over the conventional KKR-CPA, in Sec. IIIII 
we present results (configurationally-averaged density of 
states) for the application of the formalism to a one- 
dimensional model. However we emphasise that the for- 
malism presented is fully tractable for realistic three- 
dimensional systems. 



II. FORMALISM 

A. The KKR-CPA 

For the sake of clarity, we begin by briefly summarising 
the idea of the conventional KKR-CPAiSii 9 . method. The 
path operator equation describing the scattering of an 
electron in a general array of non-overlapping potentials, 
centered on site positions {Ri}, is given by 

T ij = t l 5 tJ +J2t l G{R zk )T k: > (1) 

where the underscore denotes a matrix in angular mo- 
mentum space with the usual angular momentum indices 
l,m and R 1 - 7 — Hi — R, is the position vector connect- 
ing sites i and j. For spherically symmetric scatterers the 
individual t-matrices are diagonal i.e. t l = t\ m 5ij8i m .i' m : 



but the structure constants G; m! ;' m '(R lfc ) are not. We 
shall be interested in situations where the lattice sites la- 
beled by i and j form a regular array, an infinite lattice, 
but the single-site t-matrices vary from site to site in a 
random fashion. If we consider the solution of Eq. (1) for 
all possible disorder configurations and then take its av- 
erage, we arrive at an averaged path operator r lJ which 
is translationally-invariant. In practice it is not compu- 
tationally feasible to average over every possible disorder 
configuration nor is an equation which determines the 
averaged path operator r u readily available in the form 
of Eq. 1. Under these circumstances a way forward is 
to follow the strategy of 'effective medium theories' JL2!i 
In the present context a useful 'effective medium' is that 
provided by an ordered array of effective scatterers de- 
scribed by the same t-matrix t . In the KKR-CPAfi2iiS 
the scattering amplitude describing these effective scat- 
terers is determined using the self-consistency condition 
that excess scattering off a single-site impurity embedded 
in such a medium be zero on the average. As mentioned 
in the introduction, whilst very successful in many appli- 
cations, being a single-site approximation the CPA fails 
to take into account the effect of non-local potential cor- 
relations due to the disorder in the environment of each 
site. 

B. Inclusion of non-local potential correlations 

The first step in going beyond the KKR-CPA is to 
define what we will call the non-local CPA (NLCPA) ef- 
fective medium by the following equation: 

? ij =tSij + ^2tG(R lk )? k3 (2) 

where a circumflex symbol denotes an NLCPA effective 
medium quantity. Here we have defined NLCPA effective 
local t-matrices t and a new effective propagator by 

G(R ij ) = G(R ij ) + a ij (3) 

This is composed of the usual free-space KKR structure 
constants G(R y ) which account for the lattice structure 
plus a translationally-invariant effective disorder term 
a y (= a\ 3 mVm ,). The matrix S y takes into account, 
in an averaged manner, the non-local correlations due to 
the disorder configurations. 

Clearly a cluster generalization of the KKR-CPA 
would involve determining the above NLCPA effective 
medium by requiring excess scattering off an impu- 
rity cluster of real potentials embedded in the effective 
medium to be zero on the average, and indeed this is the 
strategy followed in the next section. However the fact 
that q_ lJ is a translationally-invariant effective medium 
quantity means that it is diagonal in k , reciprocal space. 
In the present problem it is important to treat the theory 
in both real lattice and reciprocal lattice space consis- 
tently and crucially we shall deal with this by a coarse- 
graining procedure. 
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C. Cluster generalization of the KKR-CPA 

The first step is to consider a cluster of sites in the 
NLCPA effective medium. Denoting the sites within the 
cluster by capital letters, it can be shown (see Appendix 
A) that Eq. (2) can be rewritten as a cluster equation 
given by 
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with the effective cluster t-matrix defined as 



IK\7 KJ 



(4) 



(5) 



Eq. (4) is simply a re-arrangement of Eq. (2) so that 
the site matrix elements of all matrices involve the clus- 
ter sites only. The cluster t-matrix i c/ IJ describes the 
scattering within a cluster and r IJ takes account of all 
scatterings outside of the cluster via the effective clus- 
ter renormalised interactov^^ A IJ . Nevertheless, since 
G(R /J ) = G(R /J ) + a IJ , it is clear that the cluster t- 
matrix also includes non-local correlations between the 
cluster sites. ^ 

It should be stressed that A 1,1 describes all paths from 
the cluster site / to the cluster site J which only involve 
intermediate sites outside of the cluster. It is indepen- 
dent of the contents of the cluster and can be viewed 
as describing the effective medium from which the clus- 
ter has been removed i.e. replaced by a cavity. We may 
now define an 'impurity cluster' embedded in the effective 
medium simply by filling up this cavity with a particular 
configuration of site potentials. Clearly, the path opera- 
tor for sites /, J belonging to such an impurity cluster is 
given by 



—imp —cl,imp 



/ j —cl.imp — —imp 
K.L 



with the impurity cluster t-matrix defined by 



tlx, imp — t ! Sij 



Y / t I G(R IK )t_ KJ 



cl,imp 



(6) 



(7) 



Kjtl 



For a cluster containing N c sites each scattering accord- 
ing to t A or t B , there are 2^"= possible impurity cluster 
configurations. 

We are now in a position to generalize the usual CPA 
self-consistency condition to determine our approxima- 
tion to the exact configurationally-averaged medium de- 
scribed by t 1 J , namely the NLCPA effective medium de- 
scribed by £ IJ . This follows by considering for each con- 
figuration the impurity cluster path operator for paths 
starting and ending on the impurity cluster sites and de- 
manding that the average over all configurations be equal 
to the path operator for the NLCPA effective medium it- 
self i.e. 



(t u ) 

\ —imp I 



.IJ 



(8) 



The important point to note is that unlike many other 
self-consistent cluster theories such as the molecular co- 
herent potential approximation^LSS (MCPA) and exten- 
sions, here r 1,1 and t correspond to an effective medium 
which is invariant under translation from site to site and 
hence does not yield spurious gaps in the band structure 
which could affect the calculation of transport properties 
etc£ This is a consequence of our definition of the non- 
local correlation terms ex 3 as translationally-invariant. 
To preserve this periodicity, during each self-consistency 
cycle we must have 



:/./ 



J- / dk(? '-aOO-fipc) 



ik(R 7 -Rj) 



(9) 



This can be seen by applying the usual lattice Fourier 
transform to Eq. (2) and then considering the cluster 
sites /, J. Clearly, to carry out the above integration 
numerically we must have a specific representation of the 
function S(k) . Following the central idea of the DCA, 
we proceed by coarse-graining S(k) over the Brillouin 
zone. Unlike many other attempts, 5 this approach has 
been shown to be fully causali 



D. Implementation - coarse-graining of the 
Brillouin zone 

Due to the finite size N c of the cluster in real-space, 
the differences in distance between the cluster sites corre- 
spond to a set of N c cluster momenta {K„} in reciprocal 
space according to the relation 



N c ^ 



(10) 



These {K n } are at the centers of a set of iV c reciprocal- 
space patches which coarse-grain the first Brillouin zone 
of the lattice. Earlier it was noted that the effective clus- 
ter t-matrix also includes non-local correlation terms a IJ 
acting between the cluster sites. We use the above rela- 
tion to coarse-grain 2(k) as follows: 



C K„ 



S(K„) 



E 



S /J e -iK„(R 7 -R J ) 



(11) 



(12) 



Since the NLCPA maps the effective lattice problem 
to that of a self-consistently embedded impurity cluster 
problem, the fundamental assumption we make is that 
provided we are aiming to reproduce these correlations of 
finite range in the effective lattice, we may take the above 
coarse-grained values a(K„) to be a good approximation 
to S(k) for the effective lattice. In short, S(k) = S(K„) 
if k is in the n th reciprocal-space patch. This is discussed 
more formally in Appendix B. 
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The next step is to define the 'coarse-grained averaged' 
reciprocal-space matrix elements for the effective medium 
path operator by 

r(K„) = f dk (I - 1 - S(K„) - G(k)) _1 (13) 

where each integral is over the reciprocal-space patch 
f^K^ of volume N c /Qbz surrounding the point K„. This 
is straightforward since each S(K„) is constant within 
its coarse-graining patch. Defining transformations re- 
lating the real-space and coarse-grained reciprocal-space 

I 



which we use to replace Eq. (9), and we can now iterate 
to self-consistency until Eq. (8) and Eq. (16) are satisfied. 

Finally, note that the NLCPA reduces to the CPA for 
N c = 1. In this limit the non-local correlation terms a_ IJ 
vanish and the NLCPA effective t-matrix becomes equal 
to the usual CPA effective t-matrix. The NLCPA be- 
comes exact as N c — > oo where K„ — > k and non-local 
correlations over all length scales are treated. For clar- 
ity the full KKR-NLCPA algorithm is now summarised 
below. 



E. KKR-NLCPA Algorithm 

All real-space matrices in the algorithm are super- 
matrices (denoted by double underscores) in cluster-site 
and angular momentum space. For a particular energy 
E, 

1. Make a guess for the effective cluster t-matrix i c; IJ 
for the first iteration. Do this by placing an aver- 
age t-matrix (ATA), t = P(A)t A + P(B)t B (where 
P(A) is the probability of a site being occupied by 
an A atom), on each cluster site, and for the site to 
site propagation terms in the cluster use the free- 
space structure constants i.e. set a IJ = . 

2. Calculate a IJ using Eq. (5) i.e. 

where t 1 is the diagonal part of t ~ 1 . For the 
first iteration a /J will of course be zero. 

3. Convert the matrix elements a IJ to coarse- 
grained reciprocal space using Eq. (12). 



matrix elements by 

£ /J = ^E^ K ») elK " (Rl ~ R,) ( 14 ) 



r(K„)=]Tr /J e- 4K "( R '- Rj) (15) 
J 

means that the real-space effective medium path operator 
for sites /, J within the cluster is now given by 



(16) 



4. Calculate the coarse-grained matrix elements 
r(K„) using Eq. (13) and convert them to real 
space using Eq. (14). 

5. Calculate A /J by solving Eq. (4) i.e. 

6. Calculate z imp IJ for each impurity cluster config- 
uration using Eq. (6) and average over all 2 Nc con- 
figurations to obtain a new effective path operator 
at the cluster sites r IJ . 

7. Calculate the new cluster t-matrix t cl J by solving 
Eq. (4) using r /J above and A IJ from step 5 i.e. 

L = (r 1 +1)" 1 

8. Compare the new cluster t-matrix elements t cl IJ 
with those in step 1. If they are not equal to within 
the desired accuracy, repeat as necessary steps 2 — > 
8 using the new cluster t-matrix until convergence 
within the desired accuracy is achieved. 

Note that the integrations over the reciprocal space 
patches in step 4 only involve the inversion of a matrix 
in angular momentum space and therefore computational 
time is not significantly increased over the conventional 
KKR-CPA regardless of cluster size. This is in contrast 
to many other attempts such as the MCPA where the 
integration over the Brillouin zone requires the inversion 
of a super-matrix in cluster-site and angular momentum 
space for each value of k. We also draw attention to the 
recent work of Maier and Jarrell2£ where the DCA algo- 
rithm has been shown to converge more quickly than the 
MCPA algorithm, with corrections of order 1/L 2 where 
L is the linear size of the cluster. 



ij 



o^-z(7 ^(l- 1 -m^-Gik))- 1 ^^-^) 



5 



F. The cluster momenta {K n } 

Since the NLCPA maps the impurity cluster problem 
to the effective lattice problem in reciprocal space, it 
is important to realize that the real-space cluster must 
have periodic boundary conditions i.e. must preserve the 
translational symmetry of the lattice. Moreover, as ex- 
plained in Ref. |3, in order to obtain suitable reciprocal 
space patches centered at the cluster momenta {K„}, we 
must select the real-space cluster sites by surrounding 
them with a 'tile' of size L D (where D is the dimension) 
which preserves the full point-group symmetry of the lat- 
tice and only clusters which satisfy this requirement are 
allowed. 

The method for finding the corresponding cluster 
momenta {K„} satisfying Eq. (10) for a simple two- 
dimensional square lattice has been described in Ref. 0. 
Here we generalize this method to the case of three- 
dimensional simple cubic (sc) , body-centered cubic (bcc) , 
and face-centered cubic (fee) lattices commonly found in 
real disordered alloys. For the trivial case of iV c — 1, the 
real-space 'tiles' with the required symmetry would sim- 
ply be Wigner-Seitz cells surrounding each lattice point. 
For larger cluster sizes we may take the 'tiles' to be sim- 
ple cubes of volume L 3 for each of these lattices. The 
smallest possible cluster sizes are given by considering 
L = a (where a is the lattice constant) and this yields 
N c = 1 for sc (trivial), N c = 2 for bcc, and N c — 4 for 
fee lattices. The next set of allowed cluster sizes is given 
by considering L — 2a and this yields N c — 8, N c = 16, 
and N c = 32 for sc, bcc and fee lattices respectively. 

The next step is to consider the 'principal tiling vec- 
tors' {ai, &2, a 3 } where Zai+TOa2+na 3 for integers I, m, n 
map out the centers of the real-space 'tiles'. Applying 
the usual reciprocal space transformations of the form 
bi = 2n (a 2 X a 3 ) ~ (ai • (a 2 x a 3 )) etc. gives us princi- 
pal 'coarse-grained' reciprocal space vectors {bi, b 2 , b 3 }. 
We take N c non-equivalent vectors (i.e. do not differ by a 
reciprocal lattice vector) of the form Ihi + mb2 + nb 3 for 
integers ^,m, n to be our set of cluster momenta {K„} 
and these will satisfy Eq. (10). The reciprocal space 
patches surrounding these cluster momenta will be simple 
cubes of equal volume (2ir/L) 3 which together will fill out 
a volume the size of the first Brillouin zone i.e. (2ir/a) 3 
for sc, 2(27r/a) 3 for bcc, and 4(27r/a) 3 for fee lattices re- 
spectively. Integrating over these patches is equivalent 
to integrating over the first Brillouin zone of the lattice 
since K„ values and parts of patches lying outside of the 
first Brillouin zone can be translated through reciprocal 
lattice vectors (of the sc, bcc, or fee lattice as appropri- 
ate) to lie within the first Brillouin zone. In Tab. [I] we 
give some examples of sets of R/ and corresponding K„ 
values obtained using the above method. 



G. Short-range order 

The principal advantage of the KKR-NLCPA over the 
conventional KKR-CPA is that it can be implemented 
for alloys in which short-range ordering or clustering is 
present. To deal with this situation one must include an 
appropriate weighting for each of the 2 Na impurity clus- 
ter configurations in step 6 of the algorithm. Note that 
such short-range order will not destroy the translational 
invariance because it will be restored by the configura- 
tional averaging. However when using some method to 
weight the configurations it is important to bear in mind 
the periodic boundary conditions imposed on the cluster 
as explained in the previous section. To illustrate the 
above feature of the theory we will describe short-range 
order with a one-dimensional example later in this paper. 



H. Calculating observables 

In order to calculate observables such as the 
configurationally-averaged densities and density of 
states, we need to know the configurationally-averaged 
Green's function calculated within the NLCPA. The ex- 
pression for the Green's function before averaging is given 
by 

G(E, r„ r$) = £ Z* L (E, r^Z^E, r>.) 

LL' 

-£)zi,(£,r,)4(£,r0*j (17) 

L 

where L{= I, m) is an angular momentum index and 
rj(r^) lies within the unit cell centered at site 
Z l L (E, rj) and J l L {E, rj) are the regular and irregular so- 
lutions^ respectively of the single-site problem at site 
i. In the next section we show calculations for the 
configurationally-averaged density of states for a one- 
dimensional model and so here we demonstrate explic- 
itly how to take the configurational average of the site- 
diagonal Greens function. It should be stressed that 
unlike in calculations based on tight-binding models, in 
the present multiple scattering theory we solve for and 
describe the site to site fluctuations of the 'orbitals' 
Z' l L (E, rj). In short we can calculate the density r point 
by r point. It is this feature of the present theory which 
requires the following careful averaging procedure. 

The first step is to consider r and r' in the neighbor- 
hood of a cluster site /. Denoting a configuration of the 
remaining cluster sites by 7, as a generalization of Ref.l23l 
we first average over the subset of possible lattice struc- 
tures that leave the potential in site I fixed: 
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TABLE I: Examples of sets of R/ and K„ values for lattices of the sc, bcc and fee type 



lattice 


iVc 


cube edge 


Rj values (units of lattice constant a) 


K n values (units of ^) 


sc 


1 


L = a 


(0,0,0) 


(0,0,0) 




8 


L — 2a 


{Rf} = { (0,0,0), (0,0,1), 
(0,1,0), (0,1,1), (1,0,0), 
(1,0,1), (1,1,0), (1,1,1) } 


{K~} = { (0,0,0), (0,0,1), 
(0,1,0), (0,1,1), (1,0,0), 
(1,0,1), (1,1,0), (1,1,1) } 


bcc 


2 


L = a 


(0,0,0), (i,i,i) 


(0,0,0), (2,0,0) 




16 


L — 2a 


{Rf c }, {Rf c + (5, 5, |)} 


{Kf}, {Kf + (2,0,0)} 


fee 


A 

4 


Li — a 


(0,0, 0), (,2,0, j), (2'2' U /' ( U '2'2- ) 


(0, 0, U), (2,U, UJ, (U, 2,U), (0, U, 2) 




32 


L = 2a 


{Rf}, {Rf + (i,0,i)}, 
{Rf + (!,!, 0)}, {Rf + (0,i,i)} 


{Kf}, {Kf + (2,0,0)} 
{Kf + (0,2,0)}, {Kf + (0,0,2)} 



(G( J B,r,r'))7 = ^Z£( J B,r / ) ( £ P( 7 |J)<7&,)j >7 ) Z£,(£, r' 7 ) 

-J]^(£ ) r;)Ji(£,r' / ) (18) 

L 

I 



where (t££/)/ j7 is the path operator for paths starting 
and ending at site / conditionally averaged so that the 
potential on site I is known (to either be of type A or 
B) and the configuration of the remaining sites in the 



cluster is known to be 7 . P(j\I) is the probability that 
the configuration 7 of the remaining cluster sites occurs 
given the type of potential at site /. The final step is to 
average over the possible occupants of site / itself: 



<G(£,r,r')}/=£ 



LL' 



+P(B)ZE(E,rj) ^P(7|5)(ri£,) B , 7 ^ ZHErt) 
[P{A)Z£{E, n)j£(E, r'j) + P(B)ZE(E, r 7 ) jf (E, r' 7 )] 

L 

I 



(19) 



where P(A) and P(B) are the probablities that site I is an A atom or B atom respectively (i.e. the concentra- 



7 



tions of A and B atoms in the material). This can be rewritten as 
I 



LL' 



J2 P(A 7)Z£(E, r/ )(r£i,> A , 7 i# (E, r,) 

- 7 

+ £ 7)^f (£, rx) <rl£,) s , 7 z£, (25, r,) 

7 

■ 53 [P(^)Zf (£, r 7 ) (£7, r 7 ) + P(B)ZE(E, r f ) j£(25, rj)] 

I 



(20) 



where (A, 7) and (-B, 7) denote cluster configurations 
with an A atom and i? atom on site 7 respectively. The 
above expression is still exact at this stage, however it 
can be simplified using the NLCPA. In this approxima- 
tion, (t^£,)/ 7 is constructed using an 'impurity' cluster 
of configuration (I, 7) embedded in the NLCPA effective 
medium. By using Eq. (4) to eliminate the cluster renor- 
malised interactor from Eq. (6), this is given by 



( r LL')/,7 



+ t 



■cl^mp 



1 II 



LL' 



(21) 



where the impurity cluster t-matrix t has config- 

r J =cl,imp 

uration (J, 7) and the notation implies taking the I, / 
site and L, L' th angular momentum element of the super- 
matrix on the right hand side. 

It does not matter which cluster site is chosen to be 
site I in all the above formulae as G(E, r, r'), the result- 
ing approximation to (G(E, r, r')), is a translationally- 
invariant quantity. The density of states per site is given 
by 

p(E) = —Im G{E, r,r)dr (22) 
where the integral is over Oj, the volume of site /. 



III. RESULTS FOR ONE-DIMENSIONAL 
MODEL 

The multiple scattering theory in three dimensions en- 
visions three-dimensional potential wells of finite range 
surrounding the atomic nuclei, the famous muffin-tin po- 
tential. A very useful caricature of this realistic situa- 
tion can be constructed in one dimension. Evidently in 
this case the unit cell is a line segment and the poten- 
tial wells are described by a function V(x) of one spa- 
tial variable only. Interestingly, one may regard the sign 
of x as an angle and develop an analogue of the three- 
dimensional angular momentum expansion of the usual 
KKR theory. Since it was first formulated by Butler 24 



Aa 



o'o o o O'O 



(b) 



FIG. 1: A one-dimensional 'tile' for a one-dimensional four 
effective site (N c = 4) cluster is shown in (a). Its length L 
is 4a where a is the lattice constant. The corresponding set 
of cluster momenta (denoted by dots) and reciprocal-space 
patches of length 2n/L in relation to the first Brillouin zone 
are shown in (b). 



it has been made go od use of by a number of authors, 
for example see Refs. Emilia While it is computation- 
ally simple, as one might expect, it is formally identical to 
KKR in three dimensions. For example there are two 'an- 
gular momentum' values L — 0, 1 (and hence all NLCPA 
super-matrices in the algorithm have dimension 2 x N c ) 
and there is an explicit expression for the structu re con- 
stants. For a detailed description see Refs.nan 

For a cluster of size iV c , the one-dimensional reciprocal 
space patches are simply defined by the points 



(2n - N c )tt 



where a is the lattice constant and n = 1, . . . , N c . As 
an example, the real space 'onc-dimcnsional tiles' and 
reciprocal space patches for a four site (N c — 4) cluster 
are shown in Fig. ^ 

We have carried out extensive numerical calculations 
of the configurationally-averaged density of states for a 
one-dimensional alloy over a wide range of parameters. 
In all cases we have found that the KKR-NLCPA sys- 
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0.2 
Energy(Ry) 




0.2 
Energy (Ry) 



FIG. 2: Density of states for a one-dimensional model. Re- 
sults shown are for a pure A lattice, a pure B lattice and a 
KKR-CPA calculation for an AsqB2o alloy. 



co 
O 
Q 




0.2 
Energy(Ry) 

FIG. 3: A four site (N c = 4) supercell calculation for the 
^80^20 alloy obtained by averaging over all 2 4 possible con- 
figurations of an infinite periodic supercell containing four 
sites. 



FIG. 4: A four site (N c = 4) KKR-NLCPA calculation for 
the ^80^20 alloy. Notice the improved structure and the par- 
tial filling in of the band gaps compared with the KKR-CPA 
calculation in Fig. |5] 



co 
O 
Q 




0.2 
Energy(Ry) 

FIG. 5: An eight site (N c = 8) KKR-NLCPA calculation for 
the A80-B20 alloy. Notice the increasing density of states inside 
the CPA band gaps. 



tematically improves the density of states with increasing 
cluster size compared to the conventional KKR-CPA. 

As a simple illustration we set the lattice constant to be 
a = 6.00 a.u. and potential-well radius to be tmt = 2.25 
a.u. for each of the constituent potentials for A and B 
sites which are square potential wells of depth —1.2 Ry 
and —0.6 Ry respectively. The concentration of A sites is 
taken to be 80%. The density of states for the electrons in 
lattices of purely A sites and purely B sites together with 
the CPA result for the Agni^o alloy is shown in Fig. [21 
In Fig. |3 we show for comparison the results for a four 
site supercell calculation. This is obtained by considering 
an infinite periodic supercell containing four sites of a 
particular configuration and then averaging over all 2 
possible configurations. The supercell calculation is not 
the exact result due to the finite size of the supercell but 



it gives an indication of the type of structure to expect 
if we are to improve upon the CPA. A KKR-NLCPA 
calculation with a cluster size of four (N c = 4) is shown 
in Fig. 2| It is evident that much of the structure missing 
from the CPA calculation which can be associated with 
energy bands of particular configurations of the supercell 
is reproduced here. States also appear in the band gaps 
either side of the impurity band centered at 0.2 Ry which 
are absent in the CPA calculation. This is because the 
states near the band edges are the contributions of large 
clusters of like atoms and these cannot be dealt with 
by a single-site theory such as the CPA. To investigate 
this further, an eight site (N c — 8) calculation is shown 
in Fig. Clearly with increasing cluster size more and 
more states enter the band gaps. 

Next we illustrate the ability of the KKR-NLCPA 
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Energy (Ry) 

FIG. 6: Effects of short-range ordering (/3 = —0.1) and clus- 
tering (/3 = +0.1) on the four site (N c = 4) KKR-NLCPA 
calculation for the Aso£>20 alloy. 

to take into account the effects of short-range order. 
As mentioned earlier, the KKR-NLCPA can be imple- 
mented for arbitrary ensembles including those in which 
the occupancy of a site by an A or B atom is corre- 
lated to that of neighbouring sites. As a simple ex- 
ample we increase the probablity of an atom occupy- 
ing a site by a factor /3 if it follows a like atom and 
decrease its probability by the same factor if it follows 
an unlike atom in the one-dimensional cluster. Thus 
positive values of /3 correspond to short-range cluster- 
ing and negative values of /3 correspond to short-range 
ordering. As an example, for a four site impurity 
cluster of configuration ABBA we have P(ABBA) = 
(P{A) + (3) (P(B) - 0) {P{B) + (3) {P(A) - (3) where we 
have made use of the periodic boundary conditions im- 
posed on the cluster. In Fig. we show a four site 
(N c = 4) KKR-NLCPA calculation for the density of 
states using the same parameters as before along with 
short-range order parameter values (3 — —0.1, (3 = 0.0 
and f3 = +0.1. Peaks which increase or decrease can be 
identified with specific cluster configurations and the in- 
creases or decreases in the amplitude of the peaks are 
consistent with the increased or decreased cluster proba- 
bilities. 



with density functional theory. In the explicit calcula- 
tions the emphasis was on the improved structure in the 
density of states with increasing cluster size due to non- 
local correlations, and a simple example of the ability of 
the KKR-NLPCA to model the effects of compositional 
short-range order. 

In order for disordered systems to receive a 'first- 
principles' description, electronic density functional the- 
orjii& (DFT) needs to be combined with treatments 
of disorder. The self-consistent-field Korringa-Kohn- 
Rostocker coherent potential approximation (SCF-KKR- 
CPA) 28 - 29 approach is a 'first-pass' way at doing this and 
has been applied to a wide range of disordered alloysi 30 i 31 
It has also been adapted for the problem of itinerant mag- 
nets at finite temperatures whose 'disordered local mo- 
ment' spin fluctuations are handled using the CPAi 3 ^ 3 . 
In principle a DFT calculation of a disordered system 
should mean that separate SCF calculations are carried 
out to minimise the total energy for each disorder con- 
figuration individually and then an average taken over 
all disorder configurations. This is, of course, intractable 
and the strategy behind SCF-KKR-CPA calculations has 
been to minimise a functional for the averaged energy in 
terms of partially averaged charge (and spin) densities, 
i.e. the average of charge densities arising from all con- 
figurations which have either an A or a B atom on one 
site. From this approach, however, it is not straightfor- 
ward to include local environment effects such as elec- 
trostatic, 'charge-correlationi 3 ^ 3 ^ 3 ^ Madelung and local 
lattice displacement effects in alloys or indeed short-range 
order effects in general. For disordered alloys, the imple- 
mentation requires the solution of the Kohn-Sham equa- 
tions for electrons moving in the disordered arrays of the 
effective potentials associated with A and B sites which 
are insensitive to their environments. In the context of fi- 
nite temperature magnetism it means that the thermally 
induced spin fluctuations must be characterised as 'lo- 
cal moments' associated with single sites and not larger 
clusters. In this paper we have taken the first steps in 
showing how these serious omissions can be rectified by 
presenting a scheme within the KKR framework which 
goes systematically beyond the CPA. 

APPENDIX A 



IV. CONCLUSIONS 

We have presented the formalism for the KKR-NLCPA 
method which systematically improves upon the con- 
ventional KKR-CPA for describing disordered systems 
on the basis of a first-principles description of the crys- 
tal potential. We have demonstrated its use on a one- 
dimensional model and illustrated in detail the necessary 
coarse-graining procedure for real three-dimensional lat- 
tices. We have also shown how to calculate observable 
quantities with a view to combining the KKR-NLCPA 



Consider a finite-sized cluster C of sites in the effective 
medium. Eq. (2) for sites i, j belonging to the cluster can 
be written in the form 



t Si 



J2tG{R lk )? k i + Y,tG(W k )z k3 ( A1 ) 



where the sum over all sites k has been split into those 
involving sites k within the cluster and sites k outside of 
the cluster. It can be shown 5 - 36 that 



(A2) 



J2 G(K ik )t kj = J2 A i( £ y 
fc^c tec 
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with the effective ciuster renormalised interactor A y 
given by the locator expansion 

+ G{R ik )tG(R kl )tG(R lj ) + ■ ■ ■ (A3) 

k<£C,l<?C 

Inserting (A2) into (Al) and using the notation that clus- 
ter sites are denoted by capital letters gives 

1 IJ =tJu + E* (g(R ik ) + A ik ) t kj (A4) 

K 

which can be re-arranged as 

T IJ = t_ cl IJ + YAci IK A KL t LJ (A5) 

K,L 

to include the effective cluster t-matrix given by Eq. (5). 

APPENDIX B 

Firstly we note that the correlations between the clus- 
ter sites which we are aiming to reproduce in the effective 
lattice are of a finite range, say < L/2, where L is the 
linear size of the 'tile' surrounding the cluster sites (see 
Scc. IIIF(| . The next step is to divide or 'coarse-grain' the 
first Brillouin zone of the lattice of N sites into N c = L D 



patches of size (2ir/L) D (where D is the dimension) cen- 
tered at the cluster momenta {K„}. We now consider 
each coarse-grained value S(K„) to be the average of 
3(k) over the N/N c lattice momenta k within the patch 
surrounding the point K„: 

S(K„) = ^a(K„ + k) (Bl) 

k 

The fundamental assumption we make is to take these 
coarse-grained values S(K„) to be a good approximation 
to S(k) for the effective latticed This is because accord- 
ing to Nyquist's sampling theorem^! in order to repro- 
duce correlations of a finite range {L/2) in real space, we 
only need to sample the first Brillouin zone of the lattice 
at intervals of 2-k/L i.e. at the cluster momenta {K„}. 
The real space correlation terms S y ( 1 — <5y ) are cut off 
if the distance between i and j is outside the range of 
the cluster size. As N c — ► oo correlations over all length 
scales are treated since the maximum correlation length 
is proportional to the cluster size N c . 
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